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Abstract 

We study the behavior of several joint residence times of N independent Brownian 
particles in a disc of radius R in two dimensions. We consider: (i) the time Tjv(i) spent 
by all N particles simultaneously in the disc within the time interval [0,i]; (ii) the time 
(t) which at least m out of N particles spend together in the disc within the time 
interval [0, t]; and (iii) the time (t) which exactly m out of N particles spend together 
in the disc within the time interval [0, t]. We obtain very simple exact expressions for the 
expectations of these three residence times in the limit t — » oo. 



1 Introduction 

One of important and physically meaningful characteristics of individual Brownian particle 
trajectories is the time T(t) spent in a finite domain S within some observation time t. 
This time is referred to in the literature as the " occupation" ^ |5] or the " residence time" 
H3 lU OH El Ej . Note that contrary to the first exit or the first passage times out of 5, the 
residence time accounts for multiple exits from, and entries into S. As a matter of fact, 
the first passage times can be obtained, once the residence times are known (7j. The long 
time properties of T(t) are essentially dependent on the dimension of the embedding space. 
In one and two dimensions all moments of T(t) diverge with time, which mirrors the fact 
that here the Brownian motion is "compactly" exploring the space (5J. In two dimensions, 
in particular, the long time behavior of T(t) obeys the Kallianpur-Robbins' law which 
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states that the scaling variable T'(t), defined as T'(t) = — — , where D stands for the 

B w ' w Shit 

particle's diffusion coefficient and S is the area of the domain, is asymptotically distributed 
according to the exponential law. On the contrary, in higher dimensions, these moments tend 
to finite limiting values when t — * oo. Unfortunately, explicit expressions in these transient 
cases are not always available. A notable exception is the special case of a three-dimensional 
spherical domain, a case in which explicit values of the moments of the limiting values of the 
occupation time have been obtained [Hj. 

Recently, motivated in part by the development of new experimental techniques, such 
as, e.g., fluorescence correlation spectroscopy (FCS), which enables to register single particle 
events in many particle systems (see for example Ref.^SJ), there has been a considerable 
interest in the behavior of related properties of a set of independently diffusing particles (see 
Ref. [10; and references therein). For instance, the number of distinct sites visited in n steps 
by N random walks on a lattice the time spent together by the first j out of a total 
number N particles, before they escape from a given region [1211131 IT1| or the order statistics 
for iV diffusing particles in the trapping problem jlUj have been extensively analysed. It has 
been recognized that, despite the absence of any physical interaction between the diffusing 
particles, in many instances a highly cooperative behavior emerges. 

In this paper we study collective statistical properties of the mean residence time of N 
independent Brownian particles in a finite domain S. We focus here on the behavior in a 
two-dimensional continuum and suppose that S is a disc of radius R centered at the origin. 
We consider three kinds of occupation times: 

• The time Tjv(i) spent by all N particles simultaneously in the domain S within the time 
interval [0, t\; 

• The time (t) which at least m out of N particles spend together in the domain S 
within the time interval [0,t]; 

• The time T^it) which exactly m out of N particles spend together in the domain S 
within the time interval [0,t]. 

We obtain here very simple expressions for the mean values of these times in the limit t — > oo. 

We note that this type of functionals of random walks has not received much attention 
up to now, except for two particular cases. The first one concerns the occupation time of 
harmonically bounded Brownian particles in two dimensions |16j . In the second case the 
problem of the occupancy of a single lattice site by a concentration of random walkers has 
been analysed ^7]. Note that in both cases one finds that the moments of the residence time 
diverge as t — > oo. 
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The analysis of such functionals of Brownian trajectories might be of importance for several 
physical processes. Consider for example ./V molecules diffusing on a surface, which contains 
a receptor sensible to the presence of a given number of molecules, say m, with m < N, i.e. 
we suppose that there is a kind of sensitivity threshold. Suppose next that the activity of the 
receptor is proportional to the time during which it is active, that is to the time when there 
are at least m molecules in the vicinity of the receptor. We finally have a situation where the 
response of the receptor at time t is proportional to the time T*f defined previously. Another 
application can be found in the FCS which is, as mentioned, a single molecule microscopy 
method. In these experiments, a given region of space is illuminated by a laser beam, and 
one observes the fluorescence signal of the molecules going through that region. In order to 
observe a single molecule, it might be important to limit events corresponding to the presence 
of more than one molecule inside the illuminated region, by decreasing the extension of the 
beam The evaluation of the time spent when there are at least 2 particles out of the 

(2) 

total number of molecules inside the illuminated region, that is precisely the quantity T^ J , 
could give a quantitative measurement of these undesirable events, and thus could furnish an 
indication of the required extension of the beam. 



Consider ./V independent Brownian particles diffusing on an infinite two-dimensional plane. 
Let Tj(t) denote the positions of these particles at time i, while Dj stand for the corresponding 
diffusion coefficients, which are not necessarily equal to each other. We also assume that all 



Next, let us introduce three auxiliary indicator functions I5 (r), I5 ' ({r.,}) and Ig ({ r j}) 
which have the following properties 



2 The model. 



particles are initially at the origin, i.e. r!- = 0. 




(1) 




( 
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if at least m of N particles are in S, 
otherwise, 



(2) 



and 








I 



if exactly m of N particles are in S 



otherwise. 



(3) 
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The occupation times Tjv(i), Tft\t) and T^y (i), defined in the Introduction, can be then 
formally written as 



ft / N 

r N (t)= j Q at' \J[i s ( rj (t')) 



(4) 

Tfr\t)= I" dt' 1^ ({vjit')}) , (5) 
Jo 

and 

f^\t)=[ t dt'l^({v j (t')}) (6) 
J o 

Using the Poincare-type formulae |19j . the two last equations can be cast into a more conve- 
nient form: 

^ m) (*)=r^E(-i)^f*: 1 1 ) e (nwo)) - w 

J0 fc=m V" V l<ii<...<i fc <JV \j=l / 

f^(t) = fdf E(-i)^fM e (n 1^(^(0)1 , w 

J ° fc=m V V l<n<...<i fe <iV \j=l / 

where the sums extend over all ordered fc-uplets in {1, . . . , N}, and N > m > 2. Eqs.Q and 
rf)(t)and4 m 



show that Tf^\t) and T!y (f) are both functionals of T/v(i). Consequently, we focus our 



analysis on the behavior of (Tjv(t)). Results for (Tjy (i)) and (Tjy (£)) can be straightfor- 
wardly obtained once (T/v(t)} is known and will be presented here without derivation. 



3 Results. 

The first moment of Xjv(i), i.e. the mean time spent in 5 by iV particles simultaneously 
within the time interval [0, t], is given by 

M*) = <ZMi)> = jf* ^n(l 5 (r,(t')))^ df . (9) 

Turning to the limit i — > oo, and performing averaging over the realizations of Brownian 
motions, we obtain 

N 



UN 



roc \ N T / r2 ' 



(10) 



We take next the diffusion coefficient D\ of the first particle as the reference. Introducing 
the corresponding time scale r = R 2 
i = 1, . . . , N, we find that fijy obeys 



the corresponding time scale r = i? 2 /4Lq, as well as the dimensionless parameters Aj = -jj- 



HN = r ^ dv. (11) 

'o v(mt;) / 
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Next, expanding the product IT^Li (l — in powers of v, we arrive, by performing inte- 
grations by parts, at the following explicit expression: 

m = rf](-l) fc J2 (A n + ... + A ifc )ln(A n + ... + A lfc ). (12) 

k=l l<i 1 <i 2 <...<«fe<A r 

which does not seem to be an a priori trivial result. 

In the particular case when only two Brownian particles, N = 2, are present in the system, 
Eq, (|12|) yields the following symmetric compact expression: 



R 2 

= —r 



1 / D\\ 1 , ( D 2 
— In 1 + — + — In 1 + — 
£>i V D 2 J ^ D 2 \ Di 



(13) 



When one of the particles moves much faster than the second one, that is, for instance, when 
Dx » D 2 , Ea.lfTS]) reduces to 

»~£>-(£V (14) 



4Di \D 2 , 

Note that such a result appears to be quite reasonable from the physical point of view. To a 
first approximation, one estimates that the joint residence time is given by the joint residence 
time before the first exit time of the slow particle out of the disc. This time is of order 
t 2 = R 2 / D 2 . Meanwhile, the fast particle leaves the disc and returns back several times. The 
time spent at each return of the fast particle inside the disc is of order r = R 2 /Di, and the 
order of magnitude of the number of returns at time t 2 of the particle one inside the cercle 
is given by ln{t 2 /r) (cf. the Kallianpur-Robbins' law mentioned in the Introduction). As a 
consequence, we expect that for D\ 3> D 2 one has 

H2 ~ r In (— ) ~ — In (^pf \ , (15) 



D 1 \D 2 

which reproduces the dependence of ix 2 on the diffusion coefficients D\ and D 2 given by Eq. 

(ED- 

Consider now the special case when all of the particles have the same diffusion coefficient 
D\. Here, Eq. (|T2*)) becomes 

N 

m = rT (-!)"*( 7 )Hk), (16) 



N (} 

k=i \' 



which in the asymptotical limit N — > oo behaves as 

^ ~ hrM' (17) 

displaying a very slow decay function of N. 

As a matter of fact, this expression coincides exactly with the result of Refs. [2*Ul I21j 
obtained for the mean first exit time out of the disc by one of N particles. This is, of course, 
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not counterintuitive, since, after the departure of one particle, the probability that all N 
particles arrive altogether inside the disc goes to zero when N tends to infinity. 

Note also that the result in Ea. (|17|) is completely different from its lattice counterpart, 
e.g., from the result obtained for the mean joint occupation time of a single lattice site by 
N independent random walks. Here, for example, for continuous time random walks (with 
jump frequency uj) on a d-dimensional hypercubic lattice, one finds (see Ref . [T7j ) : 



/ Jl^iC '*'! ' 0)dt'= / e 
Jo J=1 Jo 



oo 

Nut' 



Ir 



ujt' 



Nd 



dt', (18) 



where Pj(0, t'\0,0) denotes the probability that the j walker, which starts its random walk 
at the origin, returns to the origin at time t', and Jq( z ) is the modified Bessel function. The 
large- N asymptotic behavior of hn thus follows ^S] 

» N = lL{ 1 + Wd + 4(^ + 2(^ + --) (19) 

In this case, decreases with the increase in the number of walkers N at a much faster rate 
than in continuum, Ea.(|17[). A similar result has been already reported in the context of the 
behavior of first passage times in finite systems |12j . 

We finally consider an important special case in which only one of the particles has the 
diffusion coefficient D\ , while the remaining N — 1 have the diffusion coefficient D2 - a case 
of an "impure" particle among a set of N — 1 "pure" particles. In this special case Ea. ()12j) 
becomes, explicitly, 



UN = —r 




+ ^EC-^Wfc + JV^r* "I \- ( 2 0) 



We turn now to the behavior of A% = lim (T^ (t)). Assuming that all particles have the 
same diffusion coefficient D\, we find, using Eq.([7)l. as well as the result obtained previously 
for /xjy, Eq. (|TT)|) . that 

W = \|<- 1 » <t "" > (I) < 21 > 

Changing the order of summations, we finally obtain the following result: 

/4 H = rm(^j EVi)"^ ~ ') HN ~ k). (22) 
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In particular, for m = 2 (which corresponds to the FCS example mentioned in the Introduc- 
tion), the last equation yields 

,(2) ( N 



^>=TN(N-l)ln\jj—j), (23) 

where r = R 2 /D\. The asymptotic limit N —* oo with m finite in Eq.(22) leads to 

# ~ — — r T : (24) 
m — 1 

which is compatible, for m = 2, with the result in Eq.(23) when iV — > oo. 

Lastly, for /Etfe 1 = lim (T^ a (t)), i.e. the mean limiting time spent simultaneously by 

t — >oo 

exactly m out of N particles within S, we find, using Eqs.(jSJ) and ()16|) (or the relation 

~(m) (m) (m+l)\ 

# = r £(-ir-*(iV - fc) Q ln(iV - fc), (25) 

We note parenthetically that, curiously enough, jlffi appears to be a non-monotonic function 
of m, as suggested by numerical analysis of Ea.(|25|). In the asymptotic limit N — > oo with m 
finite, Eq.(|lSJ) yields 

7 TT T ! (26) 

m(m — 1) 

which is not an apriori trivial result. 



4 Conclusions. 

To conclude, we have studied several types of joint residence times of a disc of radius R by 
N independent Brownian particles: The time Tjv(i) spent by all iV particles simultaneously 
in the domain <S within the time interval [0, i]; the time Tjy^(i) which at least m out of N 
particles spend together in the domain S within the time interval [0, t}; and finally, the time 
Tjy 1 ^ (t) which exactly m out of ./V particles spend together in the domain S within the time 
interval [0, t]. We have shown that in case when all the particles are initially located at the 
center of the disc, it is possible to obtain the mean values of such joint residence times exactly 
for arbitrary values of the particle number N. 
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